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ABSTRACT 

We study the linear magneto-hydrodynamical behaviour of a Newtonian cos- 
mology with a viscous magnetized fluid of finite conductivity and generalise the 
Jeans instability criterion. The presence of the field favors the anisotropic col- 
lapse of the fluid, which in turn leads to further magnetic amplification and to an 
enhanced current-sheet formation in the plane normal to the ambient magnetic 
field. When the currents exceed a certain threshold, the resulting electrostatic 
turbulence can dramatically amplify the resistivity of the medium (anomalous 
resistivity). This could trigger strong electric fields and subsequently the ac- 
celeration of ultra-high energy cosmic rays (UHECRs) during the formation of 
protogalactic structures. 

Subject headings: large scale structure of the universe — galaxies: formation — 
acceleration of particles — magnetic fields — turbulence 

A wide variety of astrophysical and cosmological problems are currently interpreted on 
the basis of gravitational instability. The current large-scale structure and galactic evolu- 
tion theories are some of the best known examples. Relatively few of the available studies, 
however, consider the role of magnetic fields, despite the widespread presence of the latter. 
Magnetic fields observed in galaxies and galaxy clusters are in energy equipartition with the 
gas and the cosmic rays. The origin of these fields, which can be astrophysical, cosmological 
or both, remains an unresolved issue (Kronberg 1994; Han & Wielebinski 2002). If magnetism 
has a cosmological origin, as observations of fiG fields in galaxy clusters and high redshift 
protogalaxies seem to suggest, it could have affected the evolution of universe (Grasso & 
Rubinstein 2001; Widrow 2002; Giovannini 2004). Studies of large-scale magnetic fields and 
their potential implications for the formation of the observed structure have been given by 
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several authors (see Thorn (1967); Jacobs (1968); Ruzmaikina & Ruzmakin (1971); Wasser- 
man (1978); Papadopoulos & Esposito (1982); Zeldovich et al. (1983); Adams et al (1996); 
Barrow et al. (1997); Tsagas & Barrow (1997); Jedamzik et al. (2000) for a representative, 
though incomplete, list). Most of the early treatments were Newtonian, with the relativistic 
studies making a relatively recent appearance in the literature. A common factor between 
almost all the approaches is the use of the MHD approximation, namely the assumption that 
the magnetic field is frozen into an effectively infinitely conductive cosmic medium. With 
few exceptions (Fennelly 1980; Jedamzik et al. 1998), the role of kinetic viscosity and the 
possibility of finite conductivity have been largely marginalized. Nevertheless, these aspects 
are essential for putting together a comprehensive picture of the magnetic behavior, par- 
ticularly during the nonlinear regime. In this article we consider a Newtonian expanding 
magnetized fluid and assume that both viscosity and resistivity are flnite. At flrst, we look 
into the linear evolution of small inhomogeneities in the cosmic medium and examine how 
the fleld and the fluid viscosity affect the characteristic scales of the gravitational instability. 
We then discuss the clectrodynamic properties of the collapsing fluid, the resulting mag- 
netic amplification and the formation of unstable current sheets. Central to our discussion 
is the concept of "anomalous resistivity", which is triggered by electrostatic instabilities in 
the plasma and can substantially reduce the electrical conductivity of the latter. We argue 
that such changes in the resistivity of the protogalactic medium will lead to the formation 
of strong electric fields during the galactic collapse. These fields can then accelerate the 
abundant free electrons and ions to ultra high energies. 

Let us consider an expanding, incompressible, magnetized fiuid with p — p{p), where p 
and p arc respectively the pressure and the density of the matter. This medium obeys the 
standard Newtonian MHD equations, which in comoving coordinates read 

I = -3%-'-V.(pu) (1) 

ot a a 

^ = V)-ir-^Vp+-V(/.+ - {V X B) X B + ^V'^u (2) 

at a a ap a Anap a^p 

VV = -AnGa'^p, (3) 

— * 

^ ^ -2^B + -Vx{uxB) + ^V'B, (4) 
ot a a 

V-S = 0. (5) 

In the above a is the cosmological scale factor, m is the fluid peculiar velocity (with V-t? = 0), 
Cg = dp/ dp is the square of the sound speed, is the gravitational potential, B is the magnetic 
field vector, u is the viscosity coefficient of the medium and rj is its electric resistivity. The 
system (l)-(5) accepts a homogeneous solution with p — pQ{t) oc a~^, B — Bo{t) oc and 
u — Uq — 0. This solution describes a weakly magnetized (i.e. Bq/pq <^ 1) Newtonian FRW 
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universe, which also defines our unperturbed background. 



Following Eq. (2), the magnetic effects are confined orthogonal to B (recall that [(V x 
B) X B] ■ B — 0), which ensures that there is no magnetic effect along the field's force fines. 
Given this, we align the background magnetic field along the z-axis of an orthonormal frame 
and consider its effects in the x-y plane. We do so by perturbing Eqs. (l)-(5) around the 
zero-order FRW solution so that p = Po + Pi, B = Bq + Bi, (j) = (j)o + 4>i and u y^O. Assuming 
wave-like perturbations (i.e. pi{r,t) = pi{t)e , Bi — Bi{t)e , etc) and using expressions 
(l)-(5), the time derivative of Eq. (2) gives 
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where H — a/ a is the Hubble parameter, = Bq/Ahpq is the Alfven speed squared and we 
have dropped the tildas for simplicity. Also, to reduce the algebra we have only considered 
perturbations orthogonal to the background magnetic field. The real component of the 
above provides a wave equation for the peculiar velocity vector, which describes a damped 
oscillation. In particular. 
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where the first term in the right shows the damping due to the expansion and of the fluid 
viscosity. The latter effect is scale-dependent and vanishes on large enough scales (i.e. as 
A; ^ 0). The last term in Eq. (7) demonstrates the conflict between gravity on the one hand 
and fluid pressure and viscosity on the other. On large scales gravity always wins and the 
perturbations collapse. Small wavelength fluctuations, however, oscillate. 

Accordingly, the magnetic presence adds to the supporting effects of pressure and vis- 

— * 

cosity only orthogonal to Bq. This means that the flrst scales to collapse along the magnetic 
held lines are smaller than those normal to them. The two critical wavelengths are the 
associated Jeans scales 
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orthogonal and parallel to Bq respectively. Overall, the magnetic presence induces a degree 
of anisotropy in the collapse. Note that for a pressureless, dust-like medium the Jeans length 
along Bq depends entirely on the viscosity and the Hubble rate. 
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As the collapse proceeds, one expects the gradual formation of turbulent motions within 
the magnetized medium. The associated eddy viscosity is proportional to Uturb ~ PoUilmix 
where Ui is the velocity perturbation and Imix the turbulent mixing length (e.g. see Biskamp 
(2003)). Assuming that Ui reaches values close to Ca, that the mixing length is a fraction 
of the magnetically induced Jeans length (i.e. Imix -C A_l) and given the low thermal tem- 
perature of the post-recombination universe, the characteristic scaling on the velocities is 
Cg < vH/pQ < c\. Note that we have also adopted the typical values ol B < 10~^ G, 
Po > 10"^'^ gr • cm-3 and H = lOO/i km • sec"^ • Mpc~\ where 0.4 < /i < 1. Then, 
A_l/A|| ~ V (CaPo) / (i^turbH) » 1, with A_l of the order of a (comoving) Mpc. Following 
the standard structure formation scenarios, the initial collapse of this very large structure 
will be followed by successive fragmentation into smaller scale formations with characteristic 
lengths A <^ A_l. Moreover, as we will outline next, the anisotropy of the collapse will further 
increase the magnetic field trapped into the gravitating medium. 

In the case of an almost spherically symmetric collapse, linear inhomogeneities in the 
magnetic energy density amplify in tune with those in the density of the matter so that 
SB^ oc Sp, where 6B^ = B^/Bq and 6p = pi/po (Tsagas & Barrow 1998; Tsagas & Maartens 
2000). Therefore, even within spherical symmetry, the formation of matter condensations 
in the post-recombination universe also signals the amplification of any magnetic field that 
happens to be present at the time. We have seen, however, that the generically anisotropic 
nature of the field will inevitably induce some degree of anisotropy to the collapse. Moreover, 
the magnetically induced anisotropy in the collapse will backreact and affect the evolution 
of the field itself. The magnetic evolution during the nonlinear regime of a generic, non- 
spherical protogalactic collapse, has been considered by a number of authors (Zeldovich 
1970; Zeldovich et al. 1983; Bruni ct al. 2003; Siemieniec-Ozieblo & Golda 2004; Dolag ct 
al. 1999, 2002; Roettiger et al. 1999). The approaches arc both analytical and numerical 
and agree that shearing effects increase the strength of the final field, while confining it 
to the protogalactic plane. Compared to the magnetic strengths of the spherical collapse 
scenario, the anisotropic increase of B is stronger by at least one order of magnitude. Thus, 
protogalactic structures can be endowed with magnetic fields stronger than those previously 
anticipated. 

So far we have seen how the magnetic presence modifies the way gravitational collapse 
proceeds, by changing the overall stability of the magnetized fiuid. This in turn affects the 
evolution of the field itself and can trigger a chain of nonlinear effects on certain scales. 
Next we will argue that this selective amplification of certain perturbative modes can play a 
important role during the nonlinear stages of protogalactic collapse, helping the instability 
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to reach its saturation point. The current induced by the total field B is 




(9) 



where V x B — aB and a measures the magnetic torsion (see Parker (1993) for details). 



Initially a is small. However, the subsequent fragmentation of the protogalactic cloud will 

increase V x B and strengthen the induced current. For example, a large-scale magnetic 
field with magnitude B ~ 10^^ G at the time of the collapse will lead to J ~ lO^a. The 
latter can reach appreciable strengths for reasonable values of a. 



cm/sec, nip is the ion mass and e is the electron charge, the excitation of low frequency 
electrostatic turbulence will increase the resistivity of the medium by several orders of mag- 
nitude (Galeev & Sagdeev 1984; Kulsrud 1998). Note that the typical critical current is 
very small in the early post-recombination universe (i.e. Jc ~ necs ~ 10~^° statamper/cm^). 
Consequently, very small values of a will lead to J > Jc, thus making the plasma electro- 
statically unstable. The effect, which is known as 'anomalous resistivity', can be explained 
through the development of current driven electrostatic instabilities in the plasma. The 
latter lead to the excitation of waves and oscillations of different kinds. The absorption of 
these waves by the ions is an additional way of transferring momentum from the electrons 
to the ions, along with the usual momentum loss from the former species to the latter. The 
average momentum loss by the electron per unit time can be written as an effective collision 
term in the form nmeVeffUe = —Ffr, where Ffr is the average friction force and n = p/nip 
the ambient number density of the plasma particles. The friction force is proportional to 
the linear growth rate of the electrostatic waves (7^) and the energy of the excited waves 
(Wjfc). The effective collision frequency is estimated to be Ueff — uJe{Wsat/kBT). Then, the 
anomalous resistivity will be 



where uje = 5.6 x 10^\/n sec~^ is the plasma frequency, ks the Boltzmann constant and 
(Wsat/kBT) ~ 1 is the saturated level of the electrostatic waves (Galeev & Sagdeev 1984). 
For certain types of current driven waves, the anomalous resistivity is several orders of mag- 
nitude above the classical one, as confirmed in numerous laboratory experiments (Bamberger 
& Friedman 1968; Yamada et al. 1975). 

This sudden switch to high electrical resistivity will inevitably lead to the formation 
of strong electric currents and therefore to a fast magnetic dissipation and intense plasma 
heating. The electric fields induced by the gravitational collapse will be £■ ~ Ca,B/c + 
VanJc ~ VanJc, givcu that Ca/c <^ 1. Thus, in this scenario, the gravitational collapse of the 



When the current exceeds the critical value Jc ~ p{e/mp)Cs, where 
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magnetized, post-recombination cloud amplifies the magnetic field and indirectly generates 
strong electric currents localized on the protogalactic plane. The anisotropy of the collapse 
enhances the local currents further and eventually drives the resistivity towards anomalously 
high values. The inevitable result is strong electric fields accelerating the abundant free 
electrons. The energy gain by an electron travelling a length A ~ A_l is Wkin ^ eEX± ~ 
GVanJc^± and the relativistic factor 7 = [1 — (w/c)^]^^ is given by 

Wkin ^Van 

JcX± e'^rjannCsX± e^By/T 
nieC^ ruec^ nieC^ mempc'^VGy/n 

in CGS units. Recall that Jc ~ encg, Cg ~ lO'^y/T and that A_l ~ B/nrripy/G when Cg < 
uH/p < cl (see (8b)). Also, we have set Wsat/ksT ~ 1 in Eq. (10), which means that 
Van ~ 1/10^ -y/n. Accordingly, the typical energy gain by a free electron can reach extremely 
high values within short timescales {t 

acc X±/c ~ 10^ yjcs^i even for relatively weak magnetic 
fields. Clearly, one can extend this process to proton acceleration and show that protogalactic 
collapse can also source UHECRs. 

We also anticipate a few particles drifting in and out these "primordial" current sheets 
(and the associated strong E- fields). If fragmentation has already taken place these par- 
ticles will diffuse along the different current sheets and possibly form the observed power 
law distribution. The detailed acceleration processes, however, is beyond the scope of this 
article (see (Azner & Vlahos 2004; Vlahos et al. 2004) for the diffusion of particles in many 
acceleration sites). 

The role of unstable currents sheets along the giant radio galaxies on the acceleration of 
cosmic ray acceleration has already been pointed out in the literature (Kronbcrg ct al. 2004; 
Colgate ct al. 2001; Nodes et al. 2003). Particles gain and lose (through synchrotron and 
inverse Compton emission) energy continuously, by travelling at speeds close to the speed 
of light. The suggestion made here is that ultra high energy cosmic-ray acceleration and 
propagation may have started almost simultaneously with the formation of galaxies through 
the electrodynamic characteristics of the gravitational instability and continue, through the 
same processes till today, since the previously described instability is active on all cosmic 
scales. It is also worth pointing out that the anomalous resistivity mechanism can easily dis- 
sipate strong galactic magnetic fields, in the form of bursty heating and particle acceleration, 
whenever rjan'^^B > V x {u x B). 

The role of cosmic magnetism during the early evolution of the first structures in our 
universe has been a subject of research and debate for many decades. Most of the available 
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studies, however, operate within the hmits of the MHD approximation, that is they assume 
a highly conducting cosmic medium. As a result, the potential large-scale implications of 
a magnetic presence within a resistive environment arc still relatively uncertain. In this 
letter, have considered a simple scenario which starts from the gravitational instability of 
a Newtonian, expanding, magnetized, viscous and resistive fluid and discusses the implica- 
tions of the field's presence during the early phases of what one might call the mild nonlinear 
protogalactic collapse. Focusing on the role of viscosity and especially on that of electrical 
resistivity, we have looked into the electrodynamic properties of the aforementioned gravi- 
tating medium and discussed issues such as current sheet formation, anomalous resistivity 
and particle acceleration on large scales. 

We begun by outlining the ways in which a magnetic presence and a finite fluid viscos- 
ity can alter the standard picture of gravitational instability. We then discussed how the 
preferential, anisotropic magnetic amplification will also increase the currents on the plane 
perpendicular to the main axis of the collapse. These gravitationally induced current sheets 
will in turn trigger electrostatic instabilities, which can then lead to anomalous resistivity 
values and subsequently to strong electric fields. We argue that the latter can be strong 
enough to accelerate the free electrons to ultra high energies. 

The infiuencc of magnetic fields on cosmic-ray propagation has been the subject of 
research in the past (see Sigl et al. (2004) and references therein) . To the best of our knowl- 
edge, however, this is the first time that a direct connection between gravitational instability 
and cosmic-ray acceleration has been suggested and discussed. We have outlined the basic 
features of this connection in a simple scenario that involves only standard Newtonian mag- 
netohydrodynamics. Given that we are in the post-recombination era and that the scales 
of interest are well within the horizon, we do not anticipate any general relativistic correc- 
tions. Special relativistic effects may need to be accounted for, but not before the particle 
velocities are an appreciable fraction of the light speed. Clearly, a detailed study of the 
acceleration mechanism proposed here should also consider nonlinear effects and the possi- 
ble implications of a varying electrical resistivity. The latter has been treated as a slowly 
changing variable, relative to the acceleration timescalc. In any case, the key requirement 
for this simple scenario to work is the presence of a magnetic field coherent on the scale of 
the collapsing protogalaxy. Our calculations argue that the required strength of this field 
is comparable to those observed in high redshift protogalaxies. If such magnetic fields are 
widespread, as current observations indicate, their amplification during the nonlinear regime 
of galaxy formation can trigger a range of nontrivial effects. In this letter we suggest that 
these effects can include the formation of strong large-scale current sheets and electric fields. 
The latter could act as driving sources for the cosmic rays observed in our universe today. 
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